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Abstract. We present a brief introduction to the relativistic kinetic theory of gases with emphasis 
on the underlying geometric and Hamiltonian structure of the theory. Our formalism starts with 
a discussion on the tangent bundle of a Lorentzian manifold of arbitrary dimension. Next, we 
' introduce the Poincare one-form on this bundle, from which the symplectic form and a volume form 

are constructed. Then, we define an appropriate Hamiltonian on the bundle which, together with the 
symplectic form yields the Liouville vector field. The corresponding flow, when projected onto 
' the base manifold, generates geodesic motion. Whenever the flow is restricted to energy surfaces 

corresponding to a negative value of the Hamiltonian, its projection describes a family of future- 
directed timelike geodesies. A collisionless gas is described by a distribution function on such 
an energy surface, satisfying the Liouville equation. Fibre integrals of the distribution function 
determine the particle current density and the stress-energy tensor We show that the stress-energy 
tensor satisfies the familiar energy conditions and that both the current and stress-energy tensor are 
divergence-free. 

Our discussion also includes the generalization to charged gases, a summary of the Einstein- 
Maxwell- Vlasov system in any dimensions, as well as a brief introduction to the general relativistic 
Boltzmann equation for a simple gas. 
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rri : INTRODUCTION 
O 

^ Kinetic theory of gases is an old subject with rich history. Its acceptance as a scien- 

tific theory with potential predictive power marked the revival of the atomic theory of 
_ nature proposed by Democritus and Leucippus in the ancient time. As early as 1738, 

rS ■ Daniel Bernoulli proposed that gases consist of a great number of molecules moving 

c3 i in all directions and the notion of pressure is a manifestation of their kinetic energy of 

motion. Gradually, the idea of Bernoulli developed further culminating in the formula- 
tion of the Maxwellian distribution of molecular velocities and the Maxwell-Boltzmann 
distribution at the end of the I9th century.' 

The arrival of the 20th century marks a new era for kinetic theory. Einstein's funda- 
mental 1905 paper on Brownian motion establishes the atomic structure of matter, and 
moreover the birth of special relativity set new challenges in re-formulating kinetic the- 
ory and its close relative thermodynamics, so that they are Poincare covariant theories. 
As early as 1911, Jiittner treated the equilibrium state of a special relativistic gas [2, 3, 4] 
while early formulations of (special) relativistic kinetic theory and relativistic thermo- 



For historical facts regarding the early development of kinetic theory see Ref. [1]. 



dynamics, successes, failures as well as early references can be found in the books by 
Pauli [5] and Tolman [6]. 

Synge in 1934 [7] (see also [8]) introduced the notion of the world lines of gas 
particles as the most fundamental ingredient for the description of a relativistic gas, 
and this idea led to the development of the modem generally covariant formulation 
of relativistic kinetic theory. Synge's idea led naturally to the notion of the invariant 
distribution function and a statistical description of a relativistic gas which is fully 
relativistic. 

The period after 1960 characterizes the modem development of relativistic kinetic 
theory based on the relativistic Boltzmann equation, and an early treatment can be found 
in the paper by Tauber and Weinberg [9]. An important contribution to the subject 
was the work by Israel [10] where conservation laws and the relativistic version of 
the H-theorem is presented and the important notion of the state of thermodynamical 
equilibrium in a gravitational field is clarified. Further it has been recognized in [10] 
that a perfect gas has a bulk viscosity, a purely relativistic effect. Gradually, with the 
discovery of the cosmic microwave background radiation, pulsars and quasars it has 
become evident that relativistic flows of matter are not any longer just mathematical 
curiosities but they are of relevance to astrophysics and early cosmology as well. These 
new discoveries led to further studies of relativistic kinetic theory and the formulation 
of the transient thermodynamics [11, 12, 13, 14, 15, 16]. The development of black hole 
physics and the realization that their interaction with the rest of the universe requires a 
fully general relativistic treatment, promoted relativistic kinetic theory to an important 
branch of relativistic astrophysics and cosmology. For an overview, see the recent book 
by Cercignani and Kremer [17]. 

The formulation of the Cosmic Censorship Hypothesis (CSH) led to the search of 
matter models that go beyond the traditional fluids or magneto-fluids, and here studies 
of the Einstein-Liouville, Einstein-Maxwell- Vlasov and Einstein-Boltzmann equations 
are very relevant and on the frontier of studies in mathematical relativity, see for instance 
Refs. [18, 19, 20, 21] and Ref. [22] for a recent review. 

The strong cosmic censorship hypothesis affirms that the maximal Cauchy develop- 
ment of generic initial data for the Einsteins equations should be an inextendible space- 
time. Since this development is the largest region of the spacetime which is uniquely 
determined by the initial data, CSH affirms that the time evolution of a spacetime can 
be generically fixed by giving initial data. In any study of the CSH the choice of the 
matter model is very important. For instance, dust and perfect fluid models may lead to 
the formation of shell crossing singularities and shocks even in the absence of gravity. 
Their formation obscures a clear understanding of the global spacetime structure associ- 
ated to the gravitational singularities. Kinetic theory offers a good candidate for a matter 
model that avoids these problems. Studies of solutions of the Einstein-Liouville system 
are characterized by a number of encouraging properties, see for instance [23] and [22]. 

Motivated by the above considerations, in this work we present a modem introduction 
to relativistic kinetic theory. Our work is mainly based on Synge's ideas and early work 
by Ehlers [24], however, in contrast to their work, we derive the relevant ingredients 
of the theory using a Hamiltonian formulation. As will become evident further ahead, 
our approach exhibits transparently the basic ingredients of the theory and leads to 
generalizations. 



In this work, {M,g) denotes a C°°-differentiable Lorentzian manifold of dimension 
n = l+ d, with the signature convention (—,+,+,...,+) for the metric g. Greek indices 
jU, V, a, . . . run from to J while Latin indices z, 7, k,. .. run from 1 to d, and we use the 
Einstein summation convention. ^{N) and I^{N) denote the class of C°° vector and 
fc-form fields on a differentiable manifold A'^, while ix and £x denote the interior product 
and the Lie derivative with respect to the vector field X. 



THE TANGENT BUNDLE 



In this section we begin with the definition of the tangent bundle and summarize some 
of its most important properties. Let T^M denote the vector space of all tangent vectors 
p at some event x E M. The tangent bundle of M is defined as 

TM := {{x,p) ■.xeM,pe T^M}, 

with the associated projection map n : TM — t- M, (jc, p) t-^ x. The fibre at x G M is the 
space 7t^^{x) = (x, T^M), and thus, it is isomorphic to TxM. 

Lemma 1. TM is an orientable, 2n- dimensional C°° -differentiable manifold. 

Remark: Notice that TM is orientable regardless whether M is oriented or not. 

Proof. The proof is based on the observation that a local chart {U, ^) of M defines in a 
natural way a local chart (V, i//) as follows. Let V := n^^{U) and define 



(l>{7^{x,p)), dxl{p) ,dxl{p),..., dxi (p) ) . 



{x,p) ^ (x^,x\...,x'',p^,p\...,p 

By taking an atlas (Ua.^a) of M, the corresponding local charts (V^, Va) cover TM. 
Furthermore, one can verify that the transition functions are C°°-differentiable and that 
their Jacobian matrix have positive determinant, yielding an oriented atlas of TM. □ 
We call the local coordinates ( ,p ,p \. . . ,p ) adapted local coordinates, 

^ j> and ^dx^^^ pY'"' ^^f^ p) } ^^'^ corresponding basis of 

the tangent and cotangent spaces of TM at (x, p) . Any tangent vector L E T(^x p^{TM) can 
then be expanded as 



and 



d 



d 

dp^ 



dx^ 



{x,p) 



^'=K,P)^Ll P^ = dpl^^iL). 



Likewise, a cotangent vector ft) G T^_^^-^{TM) can be expanded as 



0} = a^dx^'\(^^p^+l5^dp^\^^ 



{^,p) ' 



"^ = ^> dx^ 



ft) 



i^^p), 



dpi' 



i^,P)y 



The projection map 7t : TM — M induces a projection 7t^(x,p) '■ ^(x,p)(^^) ~^ 
through the push-forward of n, defined as {L) [g] '■= ^[^ ° ^] for a tangent vector L 
in T(^j^p^{TM) and a function g : M —^M. which is differentiable at x. It is a simple matter 
to verify that 



dx^ 



0, 



(1) 



and thus, the projection of an arbitrary vector field L e ^{TM) on TM is given by 



7r*(x,p)(^(x,p))=^^(^,p) ^ 



X^{x,p) 



dx^ 



+ P^{x,p)^ 



in adapted local coordinates. 

For the following, we consider a differentiable curve y : I ^ M,X (-)■ 7(A) on M. It 
induces a parameter-dependent //// which is defined in the following way: 

d 



Since TZoy = yii follows immediately that the tangent vectors X and X of 7 and 7 are 
related to each other by ii^iX) =X. In adapted local coordinates {x^,p^) we can expand 



X. 



r(A) 



■ X^{X) 



dx^ 



X 



y(A) 



:X^(A) 



y(A) 



where the coefficients are given by 

X^(A)=i^(A), P^{X)=x^{X), 

where x^{X) parametrizes the curve 7 in the local chart (^7, (^) of the base manifold, and 
a dot denotes differentiation with respect to X . 

As a particular example of this lift, consider the trajectory 7 : / — M of a particle of 
mass m> and charge q in an external electromagnetic field F E A^(M) . The equations 
of motion are 

VpP = qF{p), p:=A-Y(X), (2) 

where F : jr(M) ^ jr(M) is defined by g(X,F(y)) = F{X,Y) for all X,y e jr(M), 
and A is an affine parameter, normalized^ such that g(p, p) = — m^. Consider the tangent 
vector L to the associated lift y:I—^ TM. Since in adapted local coordinates x^ — p^ 
and p^ = —T^apP'^P^ +qF^vp'^, we find 



dxi^ 



+ 



[^,p) 



qF%{x)p'-r^,p{x)p'^pP 



(3) 



Notice that Eq. (2) implies that g{p,p) is constant along the trajectories, due to the antisymmetry of F. 



and this L defines a vector field on 7. By extending this construction to arbitrary curves 
one obtains a vector field L on TM, called the Liouville vector field. In the next section, 
we shall provide an alternative definition of the Liouville vector field based on Hamilto- 
nian mechanics. 



HAMILTONIAN DYNAMICS ON THE TANGENT BUNDLE 

The purpose of this section is twofold. At first, we introduce a symplectic structure 
on the tangent bundle which in turn is the backbone of our formulation of kinetic 
theory. Secondly, we introduce a natural Hamiltonian function on TM whose associated 
Hamiltonian vector field L coincides with the Liouville vector field defined in Eq. (3). 
Moreover, the symplectic structure defines a natural volume form on TM which will be 
helpful to set up integration on TM. 

In order to define the symplectic structure, we note that the spacetime metric g induces 
a natural one-form G A^{TM) on the tangent bundle, called the Poincare or the 
Lioville one-form. It is defined as 

0(x,p)(^) :=fo(p,;r,(,,p)(X)), X e T^,,p){TM), (4) 

at an arbitrary point {x,p) G TM. In terms of adapted local coordinates (x^^p^^) we may 
expand 



X=X^ 



and obtain 



dp'' 



P=P^ — 



X 



dx^ 



®{x,p) = S^v{x)p^ dx" I , (5) 

which shows that is C°°-differentiable. The symplectic form Q.^ on TM is defined as 
the two-form 

a := J0, (6) 
which is closed. In adapted local coordinates we obtain from Eq. (5), 

a. = g^v {x) dp^ A Jx^ I + ^ {x)p^ dx^ A dx'' I (,,^) . (7) 

The following proposition shows that Q.s induces a natural volume form on TM, and 
thus it is non-degenerated. 

Proposition 1. The n-fold product' 

n{n-\) 



A:=^^^-— a, A a, A...A^seA^" (TM) (8) 



satisfies A(^x,p) ^far all {x,p) G TM, and thus defines a volume form on TM. 



The choice for the normalization of A will become clear later 



Proof. Using Eq. (7) we find, in adapted local coordinates, 

n(n-l) 



A 



-g^vog^ivi ■ ■ ■ gfidVjdp^ A Jjc^o A dp^' A Jx^' A ... A dp^'d A dx^^ 



1 
1 



g^voE^i vi • • • g^,vAP^ ^ dp^' ^■■■^ dp^" ^ dx"^ ^ dx"' A ... A dx""' 
g^vog^, vj . . . ^^x, v,e^^' •■•^''e''"''' A . . . A Jp^ A A . . . A 



= det(^^v)^i?P° A . . . A Jp'' A A . . . A Jjc'^, 

which is different from zero since the metric is non-degenerated. 
Finally, we introduce the Hamiltonian function 

H:TM^R,{x,p)^^g,{p,p). 



(9) 

□ 

(10) 



This Hamiltonian and the symplectic structure Qs define the associated Hamiltonian 
vector field L E ^{TM) on the tangent bundle by 



dH = a, 



(11) 



which is well-defined since Q.s is non-degenerated. In order to make contact with the 
results from the previous section we work out the components of L in an adapted local 
coordinate system (x^.p^). For this, we expand 



H{x,p) = ^g^v{x)p^p^, 



dx^ 



from which we obtain 



dH=^-^P^p'dx''+g,yP^dp\ 



(12) 



On the other hand, using Eq. (7), we find 

ii^s = guy[dp^{L)dx'' -dx''{L)dp^] + ^^p^\dx''{L)dx'' -dx^mdx"] 

ax"' 

Comparing Eqs. (12,13) we conclude = p^ and 

^^'^ ~2\dx- ^ dx^ )^ P ~ ^^'^ "'^^ ^ ' 



(13) 



and thus it follows that 



dx^ ^ dp^^- 



(14) 



In the absence of an external electromagnetic field this Hamiltonian vector field co- 
incides with the Liouville vector field on TM defined in Eq. (3). Therefore, we have 
shown: 

Theorem 1. Let y be an integral curve of the Hamiltonian vector field L defined in 
Eq. (11), and consider its projection y:=iioy onto M. Then, y is necessarily an ajfinely 
parametrized geodesies: its tangent vector, p satisfies Vpp = 0, with V denoting the 
Levi-Civita connection belonging to g. 

The ideas outlined so far can be extended to diverse physical systems. As an example, 
here we discuss the case of particles with charge q interacting with an external electro- 
magnetic and gravitational field. Remarkably, a minimal modification of the Poincare 
one-form is sufficient to generalize the previous result. The modified Poincare one-form 
is 

0(x,p)(^) := gx{p, 7t,^.,p){X)) +qA,{K,^,^p){X)), X G T(^,,p){TM), (15) 

where A E A^(M) is the electromagnetic potential. In adapted local coordinates Eq. (5) 
is replaced by 

®{x,p) = [g^iv{x)p^ +qAv{x)] dx^l^^^p) . (16) 

This modification is natural in view of the fact that in the presence of an external 
electromagnetic field the canonical momentum 11^ is related to the physical momentum 
Pn by = Pii+ qA^ . The symplectic form Q.s = d& now reads 

= g^v{x) dp^ A dx^ I + 

where F = dA is, the electromagnetic field strength. Note that this Q.^ is independent of 
the gauge choice. Furthermore, Eq. (17) shows that the volume form A defined in Eq. (8) 
is unaltered. Choosing the Hamiltonian function as in Eq. (10), the resulting Hamiltonian 
vector field L coincides with the one defined in Eq. (3). 
We end this section with a simple but useful result: 

Proposition 2. We have £lH = 0, £l^s = and f^A = 0, which implies that the 
quantities H, Q.^ and A are invariant with respect to thefiow generated by L. 

Proof. Using the Cartan identity we first find £lH = iidH = —i\fls = and £lQ.s — 
diiQ-s + itd^s = —d^H = 0. With this, £iA = follows directly from the definition in 
Eq. (8). ' □ 
Remark: Since £lA = (divA^)A it follows from this proposition that the Liouville 
vector field is divergence-free, div /^L = 0. Therefore, relative to the volume form A, the 
flow in TM generated by L is volume-preserving. 

THE MASS SHELL 

We now consider a simple gas, that is, a collection of neutral or charged, spinless 
classical particles of the same rest mass m > and the same charge q moving in a 
time-oriented background spacetime (M, g) and an external electromagnetic field F. We 



dgav 
dx^ 



.x)p- + ^F^y{x) 



dx^ A Jx^ I 



(17) 



assume that the particles interact only via binary elastic collisions idealized as a point- 
like interaction."^ Therefore, between collisions, for the uncharged case, the particles 
move along future-directed timelike geodesies of (M, g) while for the charged case they 
move along the classical trajectories determined by Eq. (2). From the tangent bundle 
point of view, the gas particles follow segments of integral curves of the Liouville vector 
field L. Since all gas particles have the same rest mass, these segments are restricted to 
a particular subset r,„ of TM referred to as the mass shell. r,„ is defined as 

Tm ■.= H ' {^-'^^ = ii^^P) e TM:lH{x,p) =g,{p,p) = -m^}, (18) 

where H : TM — )■ R is the Hamiltonian defined in Eq. (10). In this section we discuss a 
few relevant properties of the mass shell. The first important property is described in the 
following lemma. 

Lemma 2. r,„ is a {2n — \ ) -dimensional -differentiable manifold. 

Proof. We consider an arbitrary point {x,p) G T,„. Since gx{p,p) = —m^ < it follows 
that n^{L) = p 0. Consequently, 7^ and it follows from dH = Q.s{-,L) and 

the non-degeneracy of Q.s that dH(^^ p-j ^ 0. Therefore, Tm is a submanifold of TM of 
co-dimension one, and the lemma follows. □ 
For the proof of the next proposition and the definitions of the current density and 
stress-energy tensors defined in the next section, the following subset of the tangent 
space TxM at a specific event x eM is useful: 

P. := {p e T,M : g,{p,p) = -m^}. (19) 

For an arbitrary (not necessarily time-orientable) spacetime {M,g), is the union of 
two disjoint sets P^ and P^ , which may be called the "future" and the "past" mass 
hyperboloid at x, respectively. Because the spacetime metric is smooth, this distinction 
can be extended unambiguously to a small neighborhood of x. However, it can be 
extended unambiguously to the whole spacetime if and only if {M,g) is time-orientable. 
In view of the definition of P^ in Eq. (19) a useful altemative definition of F,,, is 

r,n = {{x,p):xeM,peP^}. (20) 

If {M,g) is time-oriented and connected, it follows that F,,, splits into two disjoint 
components which we refer to as "future" and "past". In fact, we have the following 
proposition. 

Proposition 3. Suppose M is connected and m> 0. Then, {M,g) is time-orientable if 
and only if F^ is disconnected, in which case it is the disjoint union of two connected 
components F+ and F^. 



For the uncharged case, the self-gravity of the gas particles will be incorporated in a self-consistent 
manner in the Einstein-Liouville system. For the charged case, the self-gravity and self-electromagnetic 
field of the gas particles will be taken into account by imposing the Einstein-Maxwell- Vlasov equations 
discussed further ahead. 



Proof. See appendix A. □ 
In Lemma 2 we proved that F,,, is a submanifold of TM. In the next lemma we show 
that the volume form A on TM defined in Eq. (8) induces a {In — l)-form Q. on the 
mass shell r,„ which is non-vanishing at every point. In particular, r,„ is oriented by the 
volume form Q.. 

Lemma 3. (i) Consider the open subset V := {(x,/?) e TM : H{x,p) < 0} C TM of 
the tangent bundle. There exists a (2n— \ )-form o onV such that for all (x^p) G V 

dH(x,p) A cr(x,p) = ^{x,p) ■ (21) 

( ii) The (2n — 1 ) -form Q on F,,,, defined by the pull-back Q. .= i*o of o with respect to 
the inclusion map I : F,,, — V C TM, is independent of the choice for o in ( i) and 
defines a volume form on F,„. 

Remark: Notice that the {2n — l)-form a is not unique since we can add to it any 
field of the form dH A /3 with a (2n — 2) -form /3. However, the pull-back of o to F,„ is 
unique. 

Proof, (i) As a first step, we show the existence of a vector field G ^(V) with 
the property that dH{N) = 1 on V. For this, consider the one-parameter group of 
diffeomorphisms q)^ :V ^V, {x,p) ^ [x, e^p), A G M, which induces the rescaling 
by the factor e^ in each fibre. Let X E ^ (V ) be the corresponding generating vector 
field. 



d ; 



A=0 



Then, it follows for each {x,p) G V that 



1=0 



Therefore, A'^ := X/{2H) yields the desired vector field. 

Next, define o := ?yvA = A(A'^, ■,■,...,■). We claim that this o satisfies Eq. (21). 
In order to verify this claim, we note that it is sufficient to show that Eq. (21) 
holds at each (x.p) G V when both sides are evaluated on a particular basis of 
r(^^)(rM). A convenient basis is constructed as follows: consider the {2n — 1)- 
dimensional submanifold H = const, through {x,p). Since A^(x,p) is transverse to 
this submanifold, it can be completed to a basis {Xi := A^(^p),X2, . . . ,X2„} of 
T(^^p^{TM) such that each Xi, z = 2, . . . , 2n, is tangent to the submanifold H = const, 
that is, dHj^^ p^ (X,) = for z = 2, . . . , 2n. The claim follows by noting that 

dH{x,p) ^ <7(x,p)(^l5^2, • • • ,^2n) = <7(x,p)(^2, ■ ■ ■ ^^In) = (^(x,p) ^ ^2, ■ ■ • ,^2«)- 

(ii) Let {x,p) G T,n C V and let the basis {Xi = N(^^p^,X2, . . . ,X2n} at {x,p) be defined 
as above. Then, 



^(x,p)(^2, • • • ,^2n) = (l*<7)(x,p)(^2, • • ■ ,^2;i) = C^(x,p)(^2, • • • ,^2«) 
= \x,p) (^1 , ^2 , ■ ■ ■ , ^2n ) , 



where we have used Eq. (21) in the last step. This shows that ^[x,p) is uniquely 
determined by A^^ ^^. Furthermore, as a consequence of Proposition 1, the right- 
hand side is different from zero which proves that ^(^x p) 7^ 0. 

□ 

For the formulation of the next result, it is important to note that the Liouville 
vector field L E J%'{TM) is tangent to r,„. This property follows from the fact that 
dH{L) — iidH = £tH — 0, see Proposition 2. Therefore, we may also regard L as a 
vector field on r„j. 

Theorem 2 (Liouville's theorem). The volume form ^ on r,„ defined in the previous 
lemma satisfies 

£L^ = {divaL)n = 0. (22) 

Proof. Let o e A2"-i(y) be as in Lemma 3(i). Taking the Lie-derivative with respect to 
L on both sides of Eq. (21) we obtain, taking into account the results from Proposition 2, 



(23) 



for all {x,p) G Fm. LetX2,X3, . . . ,X2„ be vector fields on F^, and let N E ^(V) be such 
that dH{N) = 1. Evaluating both sides of Eq. (23) on {N,X2,Xti, . . .,X2n), we obtain 

= (fLC7)(X2,X3,...,Z2„) 

= L[a{X2,X3,...,X2n)]-a{[L,X2\,X3,...,X2n)-...-a{X2,X3,...,[L,X2n]) 
= L[a{X2,X3,...,X2n)]-n{[L,X2],X3,...,X2n)-...-n{X2,X3,...,[L,X2„]) 

where we have used the properties of the Lie derivative in the second and fourth step 
and the definition of Q. plus the fact that L is tangent to F,„ in the third step. Therefore, 
£lQ. = and the lemma follows. □ 
We close this section by introducing local coordinate charts on the mass shell F„,. 
For that, let (U,^) be a local chart of {M,g) with corresponding local coordinates 



[X .X , 



such that 



dx^ 



dx^ 



dxd 



xEU, 



is timelike and 4-; 

ax' 



is a basis of T^M, with the property that for each x E U, 

1,2, ... ,d are spacelike. Let (V, Xj/) denote the local chart of TM with the corresponding 
adapted local coordinates {x^ ,p^) constructed in the proof of Lemma 1. Relative to 
these local coordinates, the mass shell is determined by 

^2 - g^,{x)p^p' = goo{x){p'f + 2gojix)p'p^+g,j{x)p'pj. 



— m 



(24) 



Therefore, the mass shell F,„ can be locally represented as those {x^,p^,p') E V^(V) C 
R^" for which p^ 



p^{x^,p') with 



Plix^p' 



gOj{x)p-' ± J [goj{x)pj]^ + [-goo W] [m^ + gij{x)p'p 



-goo{x) 



(25) 



Since goo{x) < and gij{x)p'pj > for all x eU,it follows that p\{x^,p^) is positive 
and p'^{x^,p') negative. This liberty in the choice of p^ expresses the fact that locally, 
r„, has two disconnected components, representing "future" and "past". If (M, g) is time- 
orientable, this distinction can be made globally, and in this case Pj^{x^,p^) parametrize 
locally the two disconnected components Tf^ of the mass shell, see Proposition 3. 

In terms of the local coordinates of T,„, we can evaluate the volume form 

Q. = l*(?yvA) defined in Lemma 3. For that we note that by virtue of Eq. (12) for each 
(x^p) G r,„ ny the tangent vector 



1 



^^"'^^ ■ P±o{x^,P') dp^ 



with 

P±0{X^:P') ■■= g00{x)pl{x^',p')+g0j{x)pj 



= T^J[gOj{x)pj? + [-goo{x)] [m^+gij{x)p'pj], (26) 
satisfies dH{N) = 1. Therefore, we find for all {x,p) G flV, 

SI, . = i*(iNA)(^ n] = ^^^^f^^^^Jpl A . . . A A Jjc'^ A Jjc^ A . . . A dx"^ , (27) 
where in the last step we have used the coordinate expression in Eq. (9). 

INTEGRATION OVER THE MASS SHELL 

In the previous section we have introduced the mass shell F„, and the volume form 
D. on F„;. This volume form is induced from the volume form A on the tangent bundle 
TM, which, in turn is constructed from the Poincare one-form. From now on, we assume 
{M,g) to be time-oriented, such that the mass shell splits into two components T^^, see 
Propositon 3. In the following, we restrict ourselves to the "future" component r+, a 
choice which incorporates the idea that gas particles move in the future direction. 

In this section we first discuss the integral of functions defined on the mass shell r+. 
However, for the purpose of the interpretation of kinetic theory, it is also essential to 
introduce the integral of real-valued functions defined on 2(^-dimensional submanifolds 

ofr+. 

Since on r+ is defined the volume form Q., the integral of any real-valued C°°- 
functions / : r+ — )■ M of compact support is 



/a. 



r, 



+ 



For the purpose of the following analysis, we consider the particular subsets of r+ which 
are of the form 

V:={{x,p):x^K,p^P+}, (28) 



with K C M a compact subset of M which we assume to be contained inside a local chart 
{U,^) of M. Relative to such a local chart and the induced local coordinates (jc^,p') of 
r+ introduced in the previous section, the integral of a C^-function / : r+ — M of 
compact support over V takes the form 



where in these expressions, 




det(g^v(^))^rf^ 



p+o{x^,p'] 



■dct{g^y{x))d"x (29) 



X :- 



P-- 



p'+ix^p') 



dx^ 



(30) 



Provided {M,g) is oriented, the integral over p can be interpreted as a fibre integral 
over P+. This follows from the observation that is a submanifold of TxM, and this 
linear space carries the natural volume form rjx = a/ — det (gjj^vix)) dx^ A dx^ A ... A dx^ . 
Proceeding as in Lemma 3, the volume form r\j^ induces a volume form Kx on . In 
terms of coordinates {p^-,p^) of T^M chosen such that -j^ is timelike, ^ are spacelike 
for J = 1 , . . . , J, and such that p^ >Q onP^ , takes the form 



V^-det(g^v(;c)) 1 . 
P+o{x^,P') 



(31) 



with the function p^o{x^,p^) given in Eq. (26). In particular, if the coordinates {p^ ,p^) 
are chosen such that they determine an inertial frame in T^M, then we have g^v{x) = 
r\^y, and it follows that 

_dp^ ^...^dp'^ 
^JnP- + Sijp'pi ' 

from which we recognize the special-relativistic Lorentz- invariant volume-form on the 
future mass hyperboloid. 

With the definition in Eq. (31) we can rewrite the integral over p in Eq. (29) as a fibre 
integral over P+, and thus 



A-/ / 



\ 



f{^,p)7tx 



-det{g^v{x))d"x- 



K 



\ 



(32) 



where we have used again the fact that (M, g) is oriented and the definition of the natural 
volume form rj . We summarize this important result in the following lemma. 



Lemma 4 (Local splitting I). Suppose {M^g) is oriented and time-oriented. Let V C 
be a subset of the future mass shell which is of the form given in Eq. (28), where K (ZM 
is a compact subset ofM, contained in a coordinate neighborhood. Suppose / : r+ — )■ M 
is C°° -differentiable and has compact support. Then, 



( 



\ 



K 



(33) 



J 



where Kx is the fibre volume element defined in Eq. (31) and r\ is the natural volume 
element of {M,g). 

We now consider the integration of functions / : r+ — t- R on suitable submanifolds E 
of r+. The submanifolds we are considering are 2J-dimensional and oriented. In order 
to define the integral of / on E we also need a 2(i-form. Given the volume form and 
the Liouville vector field L on r,„ a natural definition for such a 2 J-form is 



(34) 



Denoting by i : E — r+ the inclusion map, it follows that for each C°°-function / : r+ — )■ 
M with compact support the integral 



i*ifco) 



(35) 



is well-defined. Note that for the particular case that L is transverse to E at every point 
of E, it follows that the pull-back of to = iiQ. is nowhere vanishing so that it defines a 
volume form on E, and in this case E is automatically oriented. We will come back to 
this case at the end of this section. 

Relative to the local coordinates {x^,p') of r+ introduced in the previous section, we 
find 



(0 



+ 



det(gpv (x)) 
1 



~7y { ^^^^^j P^g/xvi-v,^^'^' A.-.Adx"' Adp^ A...Adp" 



qf'v (x) p^ - T'ap (x) p'^p^ £ij,,„j^dpj^ A ... A dpj^ Adx^A...Adx^ 



(d-iy. 

where we have used Eqs. (3) and (27). In terms of the quantities 

1 



{n-\)\ 



det(gpv W)eAivi...v/^''' A . . . A Jx''^ 



1 V-det(g^vW) 



eij2...j^dp^^ A... A dp 



jd 



(d-iy. \p+o{x^y 

we can rewrite the right-hand side in more compact form, 

co = p^r]^An+ Lrv(jc)p^ -r«^(x)p V 



KiAl]. 



(36) 



For the particular case of hypersurfaces E of the form 



E:={(x,p):xGS,pG/\+}, (37) 

with S C M a J-dimensional, spacelike hypersurface in M which is contained in U,^ we 
have 

/ \ 



/ i*{f(o) 



1*7?^, (38) 



V/ / 

where i : S M is the inclusion map of S in M. As a preparation of what follows 
it is worth noticing that the fibre integral yields a well-defined vector field J on the 
base manifold {M,g), see Eq. (39) below. Observing that J^rj^ = ijT] we arrive at the 
following lemma. 

Lemma 5 (Local splitting II). Suppose {M,g) is oriented and time-oriented. Let E C r+ 
be a 2d -dimensional submanifold of the mass shell which is of the form given in Eq. (37), 
where S <ZM is a d -dimensional, spacelike hypersurface ofM, contained in a coordinate 
neighborhood. Suppose / : r+ — M is C°°-differentiable and has compact support. Then, 



1 1*(/£0) = J ViijT]), := j f{x,p)p^n,, 



(39) 



where Ttx is the fibre volume element defined in Eq. (31) and r\ is the natural volume 
element of {M ,g). 

So far, we have assumed the function / to be C°°-smooth and have compact support 
on r+. Although convenient from a mathematical point of view since the requirement 
of compact support avoids convergence problems, for the physical interpretation given 
in the next section this assumption on / might be too strong. In the following we give 
a brief explanation of what we believe is the correct function space for / within the 
context of kinetic theory. 

For this, we assume spacetime {M,g) to be globally hyperbolic, so that it can be 
foliated by spacelike Cauchy surfaces St,t eM. There is a corresponding foliation of the 
future mass shell r+ given by 

Lr:={{x,p):xeSt,peP+}, t eR. (40) 

Since St is spacelike, L is everywhere transverse to E/, and as discussed above it follows 
that the pull-back of (O = iiO. to E/^ defines a volume form on E;^. Therefore, we can 
define for each real- valued, C°°-function / : E/ — t- M of compact support the integral 



j fi*co. 



(41) 



^ Notice that L is transverse to E at each point of E since S is spacelike. Therefore, E is oriented. 



In fact, the integral is also well-defined for continuous functions / : — )■ M with compact 
support. In this case, it follows from the Riesz representation theorem (see, for instance 
chapter 7 in Ref. [25]) that there exists a unique measure jj-z, on the space of Borel sets 
of Lt such that 



for all such /. With this measure at hand, one can immediately define the space 
L^(Ef,^i?/i£J. We claim that this is the appropriate function space for kinetic theory. 
Primary, it turns out that the space L^(Ef, J/i^J is, in some sense, independent of the 
Cauchy surface St. The precise statement is contained in the next lemma, whose proof 
is sketched in Appendix B. 

Proposition 4. Let Ei and E2 be Id -dimensional submanifolds ofT^ with the property 
that each integral curve ofL intersects Ei and E2 exactly at one point. Then, there exists 
a diffeomorphism 9 : Ei — )■ E2 such that 



for all C°° -functions / : E2 — ?■ M of compact support. 

Secondly, the physical interpretation that the total number of particles in the gas is 
finite is accommodated in our assumption that the integral of / over E^ is finite, as we 
will see below. 

After these remarks concerning the integration of real-valued functions on the mass 
shell we are now ready to discuss physical applications of the above mathematical 
formalism. 



DISTRIBUTION FUNCTION AND LIOUVILLE EQUATION 



As we have mentioned earlier on, for a simple gas, the world lines of the gas particles 
between collisions are described by integral curves of the Liouville vector field L on the 
mass shell r„,. This section is devoted to the definition and properties of the one-particle 
distribution function and its associated physical observables. 

Following Ehlers [24] we consider a Gibbs ensemble of the gas on a fixed spacetime 
{M,g). The central assumption of relativistic kinetic theory is that the averaged prop- 
erties of the gas are described by a one-particle distribution function. This function is 
defined as a nonnegative function / : r+ — t- M on the future mass shell such that for any 
(2(i) -dimensional, oriented hypersurface E C r+, the quantity 



gives the ensemble average of occupied trajectories that pass through E. In particular, 
if E = dV is the boundary of an open set V C r+ in r+ with (piecewise) smooth 






(42) 



boundaries, then the quantity 



N[dV] = J i*{fco) (43) 

dv 

gives the ensemble average of the net change in number of occupied trajectories due to 
collisions in V. 

As explained at the end of the last section, the distribution function / lies in the space 
of integrable functions {L.dili) with respect to a hypersurface E of r+ to which L is 
transverse. However, for mathematical convenience, we assume in the following that / is 
C°°-smooth and is compactly supported on the future mass shell r+. This automatically 
guarantees the existence of all the integrals we write below. 

As a consequence of the interpretation of the integrals (42,43), we show in this section 
that for a coUisionless gas the distribution function / must obey the Liouville equation 

£Lf = 0. (44) 

The derivation of the Liouville equation (44) is based on the following result. 

Proposition 5. The Id-form (0 — is closed. Morevover, for any real-valued C°°- 
function f on we have 

difco) = {£Lf)n. (45) 

Proof. The proposition is a consequence of Liouville's theorem (Theorem 2) and the 
Cartan identity. Using these results one finds 

(fz,/)a = ftifn) = {diL+iLd){fn) = d{fiL^) = difco). 

In particular, for / = 1 it follows that da = 0. □ 
By means of Stokes' theorem and this proposition we can rewrite the integral in 
Eq. (43) as 

N[dV] = j difco) = J i£Lf)^. (46) 

V V 

In particular, let V = Uo<f<r^f be the tubular region obtained by letting flow a i2d)- 
dimensional hypersurface Eq in r+ to which L is transverse along the integral curves 
of L. Since by definition L is tangent to the cylindrical piece, ^ := Uo<r<T '^^t, of the 
boundary, it follows that / ^ fco = 0. Therefore, we obtain from Eqs. (42,46) 

A^[Er]-MEo] = JifLf)^, 

V 

that is, the averaged number of occupied trajectories at Ej- is equal to the averaged 
number of occupied trajectories at Eq plus the change in occupation number due to 
collisions in the region V. For a coUisionless gas the integral in Eq. (46) must be zero 
for all volumes V, and in this case the distribution function / must satisfy the Liouville 
equation 

£Lf = 0. (47) 



Using Eq. (3) we find 



(f^/)(x,p) =p^^(x,7,)+ [^F%Wp^-r^^WpVl ^{x,p)=0 (48) 



in adapted local coordinates of TM. The above Liouville equation can also be rewritten 
as an equation for the function f{x^,p') := f{x,p), depending on the local coordinates 
{x^ ,p^) of the future mass shell r+ introduced earlier, where we recall the definitions of 
X and p in Eq. (30). By differentiating both sides of Eq. (24) one finds 

dp\{x^,pi) _ 1 dp\{x^.P^) _ 1 ./3 

dx^ 2p+Q{x^,pjf ^ dx^ ' dp' p+oix^.pjf''^^ ' 

where here p^ = p\{x^ ,p-^) and p^ = pK Using this and 

df _ df ^ df dpi _ df p^pP dg^p df 



dx^ dx^ dpO dx^ dxi^ 2p+o dx^ dp^' 

= ^/ ^ df dp\ _ df g,ppP df 
dp' dp' dp^ dp' dp' pj^o dp^' 

Eq. (48) can be rewritten in terms of the local coordinates (x^.p^) on the mass shell as 
p'^{x^p^) + \qF\{x)p'-r^p{x)p''p^] ^{x^,p') = 0. (49) 



J dp' 



CURRENT DENSITY AND STRESS-ENERGY TENSOR 

In this section, using the distribution function /, we construct important tensor fields 
on the spacetime manifold (M, g) by integrating suitable geometric quantities involving 
/ over the future mass hyperboloidal P+. At this point, we remind the reader that {M,g) 
is assumed to be oriented and time-oriented, and thus carries the natural volume form 
Tlx defined in Eq. (31). 

We restrict our consideration to the quantities defined by: 

Jx{oc) := I f{x,p)a{p)7tx (current density), (50) 
Tx{a,l5) := J f{x,p)a{p)j3{p)7tx (stress-energy tensor), (51) 



where here a,/3 G T*M. Since / : r+ — )■ M is assumed to be smooth and compactly 
supported, these quantities are well-defined. Moreover, by construction, Jx is linear in a 
and Tx bilinear in (a,/3), and thus J and T define a vector field and a symmetric, con- 
travariant tensor field on M, respectively. Their components relative to local coordinates 



{x^') of M are obtained from J^'{x) = J^{dx^), T^^{x) = T^{dx^ ,dx^), which yields 

J^{x) - f f{x,p)p^K,= f f{x,p)- f^^J-Aa{g^,{x))d'p, (52) 

J J \pQ+{x^,p')\ V 

P+ R'' 

= / f{x,p)p^p''TZ, = [ fix J) f J-dctig^yix))d^p, (53) 



where here (jc'^ , />' ) are the adapted local coordinates of r+ with the function po+ {^^ , ) 
defined in Eq. (26). It follows from these expressions that J and T are C°°-smooth tensor 
fields on M. 

The physical significance of these tensor fields can be understood by considering 
a future-directed observer in {M,g) whose four- velocity is u. At any event x along 
the observer's world line, we choose an orthonormal frame so that g^v{x) = ^^v and 
= 5q . Relative to this rest frame we have the usual relations from special relativity, 

(p^) = (£,;.') =m7(l,v'), 7= ^ 



1 - 5ijVvj 



with v' the three-velocity of a gas particle measured by the observer at x. Therefore, in 
the observer's rest frame, we have the following expressions which are familiar from the 
non-relativistic kinetic theory of gases: 

7°(x) = J f{x,p)d''p = n{x), (particle density) (54) 
T^(x) = J f{x,p)Ed''p = n{x) <E>jc, (energy density) (56) 

R'' 

T^-'{x) — J f{x,p)p-'d'^p = n{x) < p-' >x, (momentum density) (57) 

R"' 

r'-'(jc) = j f{x,p)v^p^d'^p = n{x) <v'p^ (kinetic pressure tensor), (58) 

R"' 

where the average < A >x at x G M of a function A : F,,, — M on phase space is defined 
as 

j A{x,p)f{x,p)d'^p 

1 f id. 



We also define the mean kinetic pressure by 

p{x) := ^T'i{x) = ^n{x) < v'pi >x, 
d d 



and write p{x) := mn{x) for the rest mass density. We stress that these quantities are 
observer-dependent; so far, only the tensor fields J and T defined in Eqs. (50,51) have a 
covariant meaning. It is worth noting the following Lemma (cf. Eq. (1 11) in Ref. [24]) 

Lemma 6. The (observer-dependent) quantities £{x) := n{x) < E >x, p{x) and p{x) 
satisfy the following inequalities: 



0<dp{x)<^p{x) + 



d 



+ p{x)^<£{x)<p{x)+dp{x). (59) 



for all xEM. 



Proof The first two inequalities are obvious since p >0. For the last two inequalities, 
we first observe that 



e{x)-dp{x)=m J Y ^f{x,p)d^p. 



Since 7 ^ < 1 the fourth inequality follows immediately. Next, using the Cauchy- 
Schwarz inequality: 

Pixf = m^(jr'/-r'^'n^,p)d''p] 

< J mYf{x,p)d''p J my-^f{x,p)d'^p = £{x) [£{x)-dp{x)] , 

from which the third inequality follows. □ 
After these estimates concerning observer-dependent quantities, we now turn to co- 
variant considerations. For this, the following lemma is useful (cf. [26]): 

Lemma 7. Consider a vector p and a covariant, symmetric tensor T on a finite- 
dimensional vector space V with Lorentz metric g. Then, the following statements hold: 

(i) p EV is future-directed timelike if and only if g{p,k) < for all non-vanishing, 
future-directed causal vectors k eV. 

( ii) Suppose that T{k, k) > Ofor all non-vanishing null vectors kEV. Then, there exists 
a timelike vector k* such that 

T{k\-) = -Xg{k\.) 

with X eM.. If, in addition T{k, k) > Ofor all vectors k eV, then the timelike vector 
k* is unique up to a rescaling. 

Proof. See Appendix C. □ 
Using this lemma, we establish that at any given event x E M with /(x, ■) not iden- 
tically zero, the current density and the stress-energy tensor T^ satisfy the follow- 
ing properties: By appealing to statement (i) of this Lemma and the definition of 



in Eq. (50) we conclude that gx{Jx^k) < for all non-vanishing, future-directed causal 
vectors k G TxM. Therefore, using again Lemma 7(i), it follows that is future-directed 
timelike. Likewise, the definition of in Eq. (51) and Lemma 7(i) imply that Tx{k^ k)>0 
for all covectors k G T*M, with strict inequality if k is non- vanishing causal. Therefore, 
by Lemma 7(ii), the linear map fx : TxM T^M associated to Tx has a unique timelike 
eigenvector k^P^, 

Tx{k) = -Xk. (60) 
As a consequence, the stress-energy tensor admits the unique decomposition 

Tx = Xu®u + n, g{u,u) = —\, 7t{u,-) = 0, (61) 

where u : = k/m. Physical observers moving with this four- velocity u measure density 
e{x) = A, vanishing momentum density, and stresses described by n. For this rea- 
son, u is referred to as the dynamical mean velocity of the gas. The stress tensor n 
is orthogonal to u and symmetric; therefore, it is diagonalizable. We call the eigen- 
values Pi (x),p2 • • -iPdi^) the principal pressures. It follows from Tx(k,k) > and 
trace(r^) < that pj{x) > for all j =\,2,...,d and that £{x) > p\ [x] + p2(^) + • • • + 
Pd{^) > which implies that Tx satisfies the weak, the strong, and the dominant energy 
conditions, see Section 9.2 in Ref. [27]. 

After having introduced the current density and the stress- energy tensor, defined 
as the first- and second momenta of the distribution function over the fibre, we ask 
ourselves whether or not they satisfy the required conservation laws, provided the 
Liouville equation (47) holds. The answer is in the affirmative: 

Proposition 6. Let f : F,^ -^M.be a -function of compact support on the future mass 
shell. Then, the following identities hold for all x G M: 

divJx = J £Lf{x,p)7rx, (62) 
divTxiP) = l{£Lf{x,p))l5{p)nx + qP{Fx{J)), (63) 

for all /3 G T*M, where L is the Liouville vector field defined in Eq. (3), and F : J^T (A/) — )■ 
^{M) was defined just after Eq. (2). 

Proof. Choose a local chart {U, 0) of M and let ^ C ?7 be a compact, oriented subset 
with C°° -boundary dK in M. Consider the subset 

V:={{x,p):xeK,peP+}cr+, 

cf. Eq. (28), whose boundary is given by the 2(i-dimensional, oriented submanifold 



dV = {{x,p):xedK,peP+} 



of r+. We compute the averaged number of collisions inside V in two different ways. 
First, using Eq. (46) and the local splitting result in Lemma 4, we have 



N{dV) = j {£Lm = I I £Lf{x,p)K, 



7]. (64) 



On the other hand, using the definition of N{dV), the local splitting result in Lemma 5, 
Stokes' theorem and Cartan's identity, we find 



J^{dV) = j fOJ= j ij7]= j dijT] = j£jT]= j (div7)77. 
dv dK K K K 



(65) 



Comparing Eqs. (64,65), and taking into account that K can be chosen arbitrarily small, 
yields the first claim, Eq. (62) of the lemma. 

For the second identity, we fix a one-form /3 on M, and replace the four-current 
density 7 by 7 := r(-,/3) in Eq. (62), which is equivalent to formally replace f{x,p) 
by /(jc, p) /3 (p) in the calculation above. Then, we obtain 

div7, = j£L[nx,p)l5{p)]n, = j {£Lf{x,p)mp)Tt, + J fix,p)[£Lil5{p))]n,. 

Pt Pt 

Using local coordinates, we find on one hand div7 = V|i(r'^^/3v) = (divr)(j8) -|- 
r^^V^jSy, and on the other hand 

£LmP))=Lm^)p^]=P^P"'^l^^v+qF\{x)p'h^ 

where we have used the coordinate expression (3) for the Liouville vector field L. These 
observations, together with the definitions in Eqs. (52,53) of 7^ and T^^ yield the 
desired result. □ 



THE EINSTEIN-MAXWELL-VLASOV SYSTEM 

In this section we consider a simple gas of charged particles and take into account 
its self-gravitation and the self-electromagnetic field generated by the charges. This 
system is described by an oriented and time-oriented spacetime manifold (M, g) with 
g the gravitational field, an electromagnetic field tensor F , and a distribution function 
/ : r+ — 7- M describing the state of the gas. These fields obey the Einstein-Maxwell- 
Vlasov equations, given by 

G^v = 8;rG;v(7;!v+7)fr), (66) 
VvF^^ = qJ^, V[^F„^]=0, (67) 
£Lf = 0, (68) 



where G denotes the Einstein tensor, Gm Newton's constant, and where 



T^'^^F^aFy^'-^g^yF^pF^P (69) 
is the stress-energy tensor associated to the electromagnetic field, 

T^v' = J f{x.p)PtiPv7^x (70) 

the stress-energy tensor associated to the gas particles, and 

7^^ = j f{x,p)p^n, (71) 

is the particles' current density. Notice that the total stress-energy tensor is divergence- 
free, as a consequence of the identity (63), Maxwell's equations (67), and the Vlasov 
equation (68). Similarly, the divergence-free character of the current density follows 
from the identity (63) and Eq. (68). For recent work related to the above system, see 
Refs. [28, 29, 30]. 

THE RELATIVISTIC BOLTZMANN EQUATION FOR A SIMPLE 

GAS 

In the previous sections we developed the theory describing a coUisionless simple 
gas. For the development of this theory, we postulated that the averaged properties of 
the gas are described by a one-particle distribution function / defined on the future mass 
shell r+ and by simple arguments concluded that this distribution function obeys the 
Liouville equation f l/ = 0. 

However, undoubtedly one of the most important equations in the kinetic theory of 
gases is the famous Boltzmann equation, and in order to complete this work, in this 
section we shall sketch the structure of this equation. Like the Liouville equation, the 
central ingredient of the Boltzmann equation is again the one-particle distribution func- 
tion / associated to the gas. However, and in sharp contrast to the previous equations, the 
Boltzmann equation describes the time evolution of a system where collisions between 
the gas particles can no longer be neglected. This occurs when the mean free path (mean 
free time) is much shorter than the characteristic length scale (time) associated with the 
system. 

In order to describe the Boltzmann equation, for simplicity we consider a simple 
charged gas i.e. a collection of spinless, classical particles which are all of the same rest 
mass m > and all have the same charge q, and which interact only via binary elastic 
collisions. For the purpose of this section we shall neglect the self-gravity and the self- 
electromagnetic field of the gas, assuming a fixed background spacetime (M, g) which is 
oriented and time-oriented. For the uncharged case, the gas particles move along future- 
directed timelike geodesies of (M, g) except at binary collisions which are idealized as 



point-like interactions. If at an event x E M a binary collision occurs, then two geodesic 
segments representing the trajectories of the gas particles end and two new ones emerge. 
The nature and properties of the emerging geodesies are described by probabilistic laws 
incorporated in the transition probability through a Lorentz scalar which is the basic 
ingredient of the collision integral. The history of the gas in the spacetime consists of a 
collection of broken future-directed geodesic segments describing the particles between 
collisions. In the charged case, the same situation occurs except that the geodesic 
segments are replaced by segments of the classical trajectories of Eq. (2). 

If PiiPi stand for the four- momenta of the incoming particles and P3,/?4 for the 
momenta of the outgoing particles at the event x, then an elastic collision obeys the 
local conservation law: 

Pl+P2= P3+P4- (72) 

In the tangent bundle description such a binary collision at x E M involves four points 
(jc,/?2), (x^Ps), ix,P4) belonging to the same mass shell r+ and additionally four 
orbits of the Liouville vector field L. In particular, the orbits through {x,pi),{x,p2) 
become unoccupied while the orbits through {x, pj),{x, P4) become occupied. This 
interchange in the occupation of the orbits of L is the effect of a binary collision as 
perceived from the mass shell r+. 

Like for the case of a coUisionless system, it is postulated that its averaged properties 
are described by a distribution function / : r+ — )• M on the mass shell r+, defined by 
Eq. (42). However, the Liouville equation is replaced by the Boltzmann equation which 
has the following form: 



flfix.p) = J J J W{p3+P4^ P + P2) 

X [f{x,P4)f{x,P3) - f{x,p)f{x,P2]K^{p4)7tx{p3)7tx{p2), (73) 

where the right-hand side is the collision integral describing the effects of binary col- 
lisions. The quantity W{p3 + p4 ^-^ p + pz) is referred to as the transition probability 
scalar, and it obeys the following symmetries: 

W{p3+p4^P + P2)=W{p4+P3^P2+P): (74) 
W{p + P2^ P3+ P4) = W{p3 +P4^P + P2)- (75) 

The first symmetry is trivial, while the second symmetry expresses microscopic re- 
versibility or, as often called, the principle of detailed balancing [10]. 
Intuitively speaking the term 

W{p3+ P4 ^ P + P2)f{X:P4)f{X:P3)7txiP4)7txiP3)^x{P2) 



PI PI Pi 



in the collision integral describes the averaged number of collisions taking place in an 
infinitesimal volume element centered around x EM whose net effect is to increase the 



averaged number of occupied orbits through {x,p). On the other hand the term 




W{p^+PA f-> p + P2)f{x,p)f{x,p2)7t^{p4)7tjc{p3)7tjc{p2) 



describes the depletion of the occupied orbits through (x, p) due to the binary scattering. 

In the previous sections we have shown that the mere existence of the distribution 
function / : r+ — )■ M leads to the construction of the current density / and stress-energy 
tensor T describing the gas, and naturally this property of / remains valid for a collision- 
dominated gas. For a coUisionless gas J and T satisfy the conservation laws div/ = 0, 
divT = qF{J), as a consequence of the Liouville equation £if = 0, see Proposition 6. 
We shall show below that these properties of / and T remain valid for the case where 
the distribution function satisfies the Boltzmann equation, Eq. (73). 

In order to prove this property, let ^ : r+ — M : {x,p) i— )■ ^'{x,p) be an arbitrary C°°- 
smooth real-valued function. Multiplying both sides of Boltzmann equation by ^ and 
integrating over yields 



PIP/P/PZ 

X [f{x,P4)f{x,P3) - f{x,p)f{x,p2]n^{pA)Ttx{pj,)nx{p2)Tt^{p), (76) 

where we have made use of the symmetry properties (74,75) for the transition probability 
scalar. In particular, the right-hand side vanishes identically if m is chosen to be a 
collision-invariant quantity. The choice ^{x,p) = 1 leads to 



which implies that J is divergence-free, see Eq. (62), while the choice ^(jc,p) = I5jc{p) 
from some /3 G (M), together with the momentum conservation law, Eq. (72), yields 



which implies that dvvT^{P^) = q^x{Fx{J)), see Eq. (63). 

However, by far the most important implication of the Boltzmann equation is the 
existence of a vector field S on {M,g) whose covariant divergence is semi-positive 
definite provided that micro-reversibility holds. In order to construct this field, we 
assume for the following that the distribution function is strictly positive, with suitable 
fall-off conditions to guarantee the convergence of the integrals below. Since we abstain 
from specifying such fall-off conditions explicitly, the following arguments should be 






taken as fomial. We choose = 1 + log(A/(x,p)), where A has been inserted 

to make the argument of the logarithm a dimensionless quantity. For this choice, the 
identity (76) yields 

l[l+\og{Af{x,p))]£Lf{x,p)n,{p) = -^l 1 1 lw{pi+p2^P3+P4) 

P:t Px P.t Px P.t 

X [l0g(/l/2) - l0g(/3/4)] [/1/2 -/3/4] 7^x{Pi)'^x{P2)'^x{P1,)'^x{Pa), 

where we have abbreviated fj := f{x,pj) for j = 1,2,3,4. The right-hand side is non- 
positive, due to the positivity ofW{pi+p2^P3+P4.) and the inequality 

{logy- log x){y-x) > 0, 

which is valid for all x,y > 0. On the other hand, noting that [1 +log{Af)]£if = 
£L[log(A/)/] we obtain, by replacing the function / with the function log(A/)/ in 
Eq. (62), 

div5, = - 1 [l+logiAfix,p))]£Lnx,p)n,ip), 
Pt 

with the entropy flux vector field 5 defined by 

Sxia) :=- j log{Af{x,p))f{x,p)a{p)K,, (77) 
P.+ 

for a G T*M. Therefore, the entropy flux vector field 5 satisfies 
divS;, = \ [ [ [ [W{pi+P2^ P3+P4) 



4 

P.^PtP^P^ 

X [l0g(/l/2) -log(/3/4)] [/1/2 -73/4] 7tx{pi)7tx{p2)7tx{p3)7tx{p4) > 0. 

The above relation expresses the famous Boltzmann H-theorem. It is beyond the scope 
of this article to provide a detailed discussion of the mathematical framework underlying 
the structure of the collision integral and related properties of the Boltzmann equation. 
The reader is referred to Refs. [31, 18, 17] for a more detailed discussion and properties 
of the Boltzmann equation. 



CONCLUSIONS 

In this work, we have presented a mathematically-oriented introduction to the field 
of relativistic kinetic theory. Our emphasis has been placed on the basic ingredients that 
constitute the foundations of the theory. As it has become clear, of a prime importance 
for the description of this theory are not any longer point particles but rather and 
according to Synge's fundamental idea, world lines of gas particles. From the tangent 
bundle point of view, these world lines appeared as the integral curves of a Hamiltonian 



vector field. It is wortli stressing here tlie liberty and the flexibility that characterizes the 
relativistic kinetic theory of gases. The Poincare one-form as well as the Hamiltonian 
and thus the structure of the Hamiltonian vector field are at our disposal. In this work 
we have made simple and natural choices for the Poincare one-form and Hamiltonian, 
and we were led from first principles to the Liouville equation, while for the case where 
the self-gravity and self-electromagnetic field of the gas are accounted for we arrived 
naturally to the Einstein-Liouville and the Einstein-Maxwell- Vlasov systems. 

Therefore, given the freedom in the Poincare one-form and the Hamiltonian it is worth 
thinking of relativistic kinetic theory models describing dark matter or dark energy. This 
could be interesting from an astrophysical and cosmological point of view. 
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APPENDIX A: PROOF OF PROPOSITION 3 

We prove that r,„ is connected if and only if {M,g) is not time-orientable. Suppose 
first that r„, is connected and let us show that {M,g) cannot be time-orientable. For 
this, choose an arbitrary point x G M and two timelike tangent vectors and 
at X such that k± E P^. Since F,,, is connected by hypothesis, there exists a curve 
7: [0, 1] — )• F,„,? {y{t),k{t)) which connects the two points {x,k^) and The 
projection of this curve yields a closed curve 7 : [0, 1] — )■ M in M through x. Along this 
curve is defined the continuous timelike vector field k{t) which connects k^ and (see 
Fig. 1). Therefore, {M,g) is not time-orientable. 

Conversely, suppose {M,g) is not time-orientable. We shall prove that F,,, is con- 
nected. In order to show this we note the hypothesis implies the following. There exists 
an event x E M,a closed curve 7 : [0, 1] — t- M through x and a continuous timelike vector 
field k{t) E Tyi^f^M along 7 which we may normalize such that g{k{t),k{t)) = —m^ for 

aU t E [0, 1], with the property that k- := k{0) E and kj^:=k{l)E P^. In this way we 
obtain a curve 7: [0, 1] r,„,t H- {Y{t),k{t)) in F,,, which connects (x^k-) with {x,k+). 
Since the sets P^ are connected, it follows that any two points (x,/?), {x,p') E F,„ in the 
same fibre over x can be connected to each other by a curve in F^. 

Now let us extend this connectedness property to arbitrary points {x2,P2) £ 

F„,. Since M is connected there exist curves 7i, 72 in M which connect x with xi and x 
with X2, respectively. By parallel transporting pi along ji and p2 along fi we obtain 
curves 71 and 72 in F^ which connect (xi.pi) with {x,p) and {x2,P2) with {x,p') 
respectively.^ By the result in the previous paragraph, {x, p) and (jc, p') can be connected 
to each other by a curve in F^. Therefore, (jci,pi) and {x2-,P2) can also be connected to 
each other and it follows that F,„ is connected. 



7i and 72 are the horizontal lifts of 71 and 72 through the points {x\,p\) and (x2,P2)> respectively. 




In order to conclude the proof we need to show that r,,, consists of two connected 
components if (M, g) is time-orientable. In order to prove this, choose x EM and a time- 
orientation P+ and at x. Given any point y G M we may connect it to x by means of 
a curve 7 in M since M is connected. We choose the time orientation P^ at y such that 
kx E P^ if and only if ky E Py for all parallel transported timelike vector fields k along 
7. This choice is independent of 7 otherwise (M, g) would not be time-orientable. In this 
way we obtain the two connected subsets 

Tt:={{x,p)ET,n:pEP^} 

of r,„ which are disjoint and whose union is r,„. 



APPENDIX B: SKETCH OF THE PROOF OF PROPOSITION 4 

We define 9 : Ei — E2 to be the map that associates to each point piE'Li the unique 
point p2 E E2 that lies on the integral curve of L through pi. By transporting the values 
of / along the integral curve of L in the segment between Ei and E2, we can extend / 
to a C°°-function / : r+ — t- M with compact support, such that £if = in the region V 
between Ei and E2. It follows by Stokes' theorem that 

j ilifco) - J ilifco) = J difco) = I {£Lm = 0, 

£2 £1 V V 

where ij : Ej — )■ r+ denote the inclusion maps for j = 1,2, and where we have used 
Proposition 5 in the second step. Since l|/ = / and l[/ = (p*f the Proposition follows. 



APPENDIX C: PROOF OF LEMMA 7 



Without loss of generality we can assume that {V,g) is Minkowski spacetime, so that 
the standard special relativistic notions of causality hold. We work in coordinates (x^) 
for which the metric has the standard diagonal form. 

(i) At first we note that any non-vanishing causal, future-directed vector k eV has 
the components {k^^) = {k^,k') with k'^ > and dijk'kJ < (k^)^. Suppose first that 
p EV is future-directed timelike. Then, after a Lorentz transformation, (p^) = 
(jc'^,0,...,0) with p^ > and it follows that g{p,k) = -p^k^ < 0. Conversely, 
suppose that g{p,k) < for all non- vanishing causal, future-directed k. After a 
rotation the components of p are (p^) — {p'^, p^ ,0, . . . ,0) and the two choices 
(fc^) = (l,l,0,...,0)and (^^) = (1,-1,0,...,0) lead to > = +pi 
and > g{p,k-) = —p^ — p^ which implies \p^ \ < p^, and hence, that p is future- 
directed timelike. 

Remark: In fact, as the proof shows, it is sufficient to require k to be non-vanishing 
null and future-directed in order to show that p is future-directed timelike. 

(ii) Consider the unit mass hyperboloid 

■={ueV : g{u,u) = -1,M^ > 0}, 
whose elements may be parametrized according to the stereographic map. 



V 



with Bi (0) := {v G M-^ : dijv^v-^ < 1} the unit open ball centered at the origin. Then, 
we have on H^, 

T{u,u) = [ToQ + lT^y + TijVv^] =: /(v), 

where the function / : 5i (0) — )■ M is smooth. Suppose v — )■ e, |e| = 1, approaches the 
boundary of its domain. Then, since Tqq + 27(),v' -|- Tijv'v-' — j- T{k,k) with the null 
vector k = (l,e), and since by assumption T{k,k) > 0, it follows that /(v) — )■ o°. 
Therefore, the function / has a global minimum at some point v* G 5i (0). Taking 
the gradient on both sides of ^ 

(1 - |v|2)/(v) = Too + lTo^ + Tijv'vj 

and evaluating at v = v*, we obtain 

-/(v*)v; = ro,-f7^-,(v*)^ 



^ Alternatively, one might also invoke the Lagrange multiplier method, applied to the restriction of the 
quadratic form r(M,M) on in order to conclude that at the minimum m*, T{u*,-) = —Xg{u*,-). 



from which we also have -/(v*)|v*p = Toi{v*y + Tij{v*y{v*y and so /(v*) = 
Too + Toi{v*y. Therefore, the timelike vector k* := (l,v*) satisfies T^v{k*y = 
-Xkl with A =/(v*). 

As for uniqueness, we first observe that after a rescaling and a Lorentz transfor- 
mation we may assume that k* = (1,0,0, ... ,0). Then, Tqq = X, Tqj = 0, and by 
applying a rotation if necessary, we can assume that Tij is diagonal, such that 
(7)xv) = diag(A,pi,/>25 • • • 5P<^)- It follows from the hypothesis that A > and 
/)j > for all 7 = 1 , 2, . . . , J. Furthermore, A = implies p j > for all j = 1 , 2, . . . , J 
because of the hypothesis. Therefore, the eigenspace belonging to the eigenvalue 
—A of the matrix (T^y) = diag{—X, pi,p2, ■ ■ ■ ^ Pd) must be one-dimensional since 
—A 7^ pj for all J = 1 , 2, . . . , J. 

Remark: Due to the fact that the scalar product defined by the Lorentz metric g is 
not positive definite, it is not always the case that the stress-energy tensor T^y can be 
diagonalized, despite of the fact that it is symmetric. In fact, the results above show that 
T^v is diagonalizable if and only if it admits a timelike eigenvector. An explicit example 
for which the condition in assumption (ii) in Lemma 7 is not satisfied is 

with a null vector k, corresponding to a stress-energy tensor for null dust. In this case, 
the only eigenvalue of T^^y is zero, and its eigenspace consists of the vectors which are 
orthogonal to k. In particular, this example shows that the imposition of the weak energy 
condition is not sufficient to guarantee the diagonalizability of the stress-energy tensor. 
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